Geant4 simulation of fast electron bremsstrahlung imaging at the HL-3 tokamak* 
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To further research on high-parameter plasma, we plan to develop a two-dimensional hard X-ray (HXR) imag- 
ing system at the HL-3 tokamak to measure HXRs with energies ranging from 20 to 300 keV. The application of 
an array-structured detector ensures that this system can measure HXR-radiation spectra from the entire plasma 
cross section. Therefore, it is suitable for the study of fast-electron physics, such as radio-frequency wave 
current drives, fast electrons driving instabilities, and plasma disruptions in fusion research. In this study, we 
develop a simulation for calculating fast-electron bremsstrahlung in the HL-3 tokamak based on the Monte Carlo 
simulation code Geant4, in which the plasma geometry and forward scattering of fast-electron bremsstrahlung 
are considered. The preliminary calculation results indicate that the HXR energy deposition on the detector 
is symmetrically distributed, even though the plasma distribution is asymmetric owing to the toroidal effect. 
These simulation results are helpful in constructing the relationship between the energy deposition on the de- 
tector and parameter distribution on the plasma cross section during HL-3 experiments. This is beneficial for 
the reconstruction of the fast-electron distribution function and for optimizing the design of the HXR-imaging 


system. 
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I. INTRODUCTION 


In tokamak plasma, fast electrons can be produced in 
many ways, such as through the application of auxiliary heat- 
ing methods, disruptive events, and sawtooth activities. Al- 
though their density is two orders of magnitude below the 
plasma density, fast electrons play an important role in main- 
taining the steady operation of tokamaks with high parame- 
ters [1, 2]. During radio-frequency (RF) wave heating and 
current-driven experiments, electrons gain energy from the 
RF waves through resonance effects and then become fast 
electrons. The redistribution of fast electrons and their colli- 
sions with other plasma ions and electrons eventually increase 
the plasma temperature and current. However, owing to their 
low Coulomb-collision resistance, they can be converted into 
runaway electrons, which form a runaway beam. Once the 
runaway beam is formed, it may affect the normal operation 
of the experiments, causing melt damage to the plasma-facing 
components and divertors [3—7]. In more serious cases, it ter- 
minates the discharges and compromises the safe operation 
of the device. Hard X-rays (HXRs) are produced from col- 
lisions between fast electrons and bulk plasma ions and pro- 
vide a considerable amount of information about fast elec- 
trons. Therefore, the measurement of HXRs during plasma 
discharge is an important method for studying fast-electron 
physics, such as the build-up of the electron-distribution func- 
tion, electron cyclotron wave-power deposition profile [1, 2], 
and dynamics of the instabilities [8, 9]. 
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Most traditional HXR diagnostics depend on arrays of de- 
tectors at different toroidal and poloidal locations around the 
tokamaks, such as HT-6B [10], HT-7 [11, 12], EAST [13], 
HL-1M [14], HL-2A [15, 16], and EXL-50 [17]. Multidi- 
rectional detection has enabled the reconstruction of emis- 
sion profiles from the plasma cross section [8]. In the TCV 
tokamak, a multichannel HXR spectroscopy system allows 
observation of the evolution of the electron distribution re- 
lated to fishbone oscillations [18]. In the WEST tokamak, 
the designed two-dimensional (2D) HXR camera provides 
spatial, temporal, and energy-resolved measurements of the 
HXR emission from the full plasma cross section, enabling 
the investigation of the electron temperature and fast-electron 
tail-density evolution, as well as the characteristics of tung- 
sten X-rays due to beam-target emission at the edge [19]. In 
the HL-2A tokamak, measurements from an HXR pinhole 
camera indicate that the spatial distribution of fast electrons 
during the lower hybrid current drive (LHCD) has a peaked 
profile, implying that fast electrons are mainly produced in 
the plasma core [20, 21]. The spatially distributed HXR de- 
tection array at the J-TEXT tokamak showed that most run- 
away electrons transported to the plasma boundary tend to be 
lost on the limiters [22]. In other fusion devices, more than 
two arrays are used for tangentially viewing HXR diagnos- 
tics, with the help of tomography algorithms; this method has 
promising applications for the analysis of MHD instabilities 
in the T-10 tokamak [23] and current-drive efficiency studies 
in Tore Supra [24]. 

Owing to the anisotropic characteristics of bremsstrahlung 
and irregular shape of plasma, inverting line-integrated mea- 
surement data to obtain a 2D emissivity distribution is dif- 
ficult [25]. Therefore, 2D HXR-imaging diagnostics have 
been developed to consider space requirements and avoid 
complex correlation and inversion procedures. In the PBX- 
M tokamak, a 2D HXR camera was applied to obtain direct 


measurements of the location of the current-driven electrons, 
which have been used for current profile-modification exper- 
iments [26-28]. Additionally, the HXR camera has the ca- 
pability to observe bremsstrahlung radiation from fast elec- 
trons during LHCD experiments, thereby allowing accurate 
fast-electron-transport studies [29, 30]. In the DIII-D toka- 
mak, an HXR camera was installed on the outer midplane of 
the tokamak with 123 detectors, which covered most of the 
plasma volume [31]. An electron-distribution function with 
2D spatial resolution was obtained during runaway-electron 
experiments, allowing the characterization of the evolution of 
the RE distribution function during RE growth and dissipa- 
tion [32-34]. 

Two-dimensional HXR measurements have been success- 
fully implemented in many tokamak devices; however, owing 
to the spatial asymmetry of the plasma induced by toroidal 
effects, associating tangential measurements with informa- 
tion on the plasma cross section is difficult. Computing the 
2D distribution function of fast electrons and studying the 
instability physics related to fast electrons is also challeng- 
ing. Therefore, tangential 2D HXR imaging systems have 
not been widely implemented in tokamaks. Under such con- 
ditions, we conducted a 2D HXR-imaging simulation study 
on the HL-3 tokamak to determine the correspondence be- 
tween the images of the tangential HXR-imaging system and 
plasma cross-section parameters. The significance of the sim- 
ulation was evident. First, the simulation results can quali- 
tatively show whether the HXR-imaging system is feasible. 
Second, the HXR-imaging system can be optimized by the 
simulation results, such as the placement position and shield- 
ing structure. Third, the simulation results can help us an- 
alyze the measured experimental results and then inversely 
deduce the radiative properties and distribution characteris- 
tics of the fast electrons. The Geant4 simulation model is 
based on the plasma geometry of the HL-3 tokamak and 
HXR imaging system, and a fast-electron forward-scattering 
model is also considered. The remainder of this paper is or- 
ganized as follows. Sec. II introduces the HL-3 tokamak and 
HXR diagnostics. A Geant4 simulation model based on the 
plasma configuration of the HL-3 tokamak and fast-electron 
bremsstrahlung model containing forward-scattering effects 
are presented in Sec. III. The preliminary simulation results 
are presented in Sec. IV. Sec. V provides potential reasons for 
the energy deposition that appears on the detector. Finally, the 
summary and discussion are presented in Sec. VI. 


Il. HL-3 TOKAMAK AND HARD X-RAY DIAGNOSTIC 


A schematic of the HL-3 tokamak is shown in Fig. 1, which 
shows the main parts of the device and the HXR-imaging sys- 
tem. The transparent purple rings denote plasma, and the yel- 
low cylinders are the sight lines of the detector arrays. The 
HL-3 tokamak is a medium-sized device with a major radius 
of R = 1.78 m, a minor radius of a = 0.65 m, an elonga- 
tion of x > 1.8, and a triangularity of ô > 0.5. The main 
purpose of this device is to support the operation of future 
fusion devices by integrating the technological and physical 


aspects [35, 36]. The designed toroidal magnetic field is ap- 
proximately 2.2-3 T, and the maximum plasma current is 3 
MA. The heating and current-drive systems can provide a to- 
tal power of 27 MW. The first plasma was successfully pro- 
duced in 2020 [37], and currently, the maximum plasma cur- 
rent is higher than 1 MA. 

The tangentially viewed HXR-imaging system primarily 
comprises an LYSO scintillator crystal, a SiPM detector ar- 
ray, and online spectra-processing electronics [38]. The sys- 
tem is located outside the vacuum chamber, and a 1-mm-thick 
beryllium window is used to protect the vacuum. This system 
can measure HXRs with energies ranging from 20 keV to 300 
keV, with a spatial resolution of several centimeters, a tem- 
poral resolution of several tens of microseconds, and an en- 
ergy resolution of approximately 8 %@662 keV. The design 
of the HXR imaging system is illustrated in Fig. 2, where the 
LYSO scintillator and SiPM array are located inside the tung- 
sten—copper shield; only photons passing through the small 
pinhole in front of the shield can enter the detector. When a 
LYSO crystal is struck by X-ray photons, certain wavelengths 
of visible or near-visible photons are scintillated. The visi- 
ble light is detected by the SiPM arrays and converted into 
electrical-pulse signals, and an analog-to-digital converter 
then digitizes the electrical-pulse signals and delivers them 
to the FPGA module. The FPGA module rearranges these 
digitized signals into a specific format to identify the photon 
energy and deposition time. Finally, the experimental data 
containing time and pulse-amplitude information are trans- 
ferred to a remote-controlled PC. The time evolution of the 
radiation intensity is obtained from the time evolution of the 
energy deposition on each detector. The spatial positions cor- 
responding to different detectors can reflect the spatial evolu- 
tion of the radiation intensity. The energy-spectrum informa- 
tion of a specific period can be used to deduce the temperature 
and distribution characteristics of fast electrons. 


Ill. INTRODUCTION OF GEANT4 SIMULATION MODEL 
AND FAST-ELECTRON BREMSSTRAHLUNG MODEL 


The primary objective of this study is to model X/y-ray 
photon-energy deposition on the detector. Similar work has 
been performed in the EAST tokamak, which simulates the 
energy response of a scintillating fiber detector using the 
Geant4 code [39, 40]. We also performed simulations in pre- 
vious studies, where fast-electron forward scattering was not 
considered [41]. Initially, we assumed that the direction of 
movement of each photon was arbitrary. After testing a large 
number of particles and a long calculation time, the deposited 
energies on the detector remained very small and did not ex- 
hibit any distributional properties. 

We performed a simple analysis to explain these results. 
Because of the shielded collimation system at the front of the 
HXR-imaging system, only incident-ray photons entering the 
collimation pinhole can deposit their energy on the detector. 
On the sight lines of the imaging system, only particles in the 
cross section of the plasma that are parallel to the surface of 
the detector have the greatest probability of interacting with 
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Fig. 1. (Color online) Three-quarter cut-away view of the HL-3 
tokamak, where the HXR-imaging system is located between the 
toroidal-field coil and vacuum chamber. The translucent purple 
shape is the plasma, and the yellow arrows are the optical path. Other 
components such as the support frame, vacuum chamber, poloidal 
field, and toroidal-field coils are also displayed in the figure. 
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Fig. 2. (Color online) Cut-away view of the HXR system. The LYSO 
scintillator and SiPM detector are placed inside the shields, and the 
optical path can cover most of the plasma cross section. 


the detector. We use this plane as an example to estimate 
the value of the maximum probability, where the distance be- 
tween the collimating pinhole of the HXR-imaging system 
and plasma cross section parallel to the detector surface is ap- 
proximately 2.5 m, and the designed collimation pinhole has 
a diameter of 2mm. We obtain the following: 


Area of the pinhole 
Area of the radiation sphere 
_ mx 0.002? l 
Axa x 2.52 


imax = 
(1) 
= 1.6 x 1077 


However, this maximum probability should be divided by 
a large number (e.g., 50) when considering the entire HL- 
3 plasma region because the farther the photon is from the 
plane parallel to fmaz, the less likely it is to be deposited 
on the detector. The probability of a photon moving in an 
arbitrary direction and depositing its energy on the detector is 
approximately: 


fmaz =0 
en =3.2 x 107°. 2 
Bp = 32x 10 (2) 


That is, for every 101° X/-+-ray photons, approximately 32 
of them will deposit their energy on the detector. To elucidate 
the energy-deposition distribution, the number of deposited 
particles would be of the order of one million; thus, the num- 
ber of X/ y-ray photon sources would be uncountable. This 
is detrimental to our calculations because this number of par- 
ticles would be impossible to enter into the Geant4 code; the 
larger the number of particles, the more resources are con- 
sumed by the calculations. 

To resolve this problem, we assume that all incident parti- 
cles move toward the collimation pinhole. Under a particle- 
source distribution with a slab configuration, the deposited 
energies on the detector are similar to the distribution of the 
electron density or temperature on the plasma cross section. 
Under a particle-source distribution with an HL-3 plasma 
configuration, the high-field side has more X /y-ray photons 
because of the longer chord length, resulting in a deviation 
of the Gaussian distribution of the energy deposition on the 
detector. As this model does not match the actual conditions, 
this part of the calculation deserves further speculation. 

In the latest studies, we updated our simulation model 
based on the HL-3 plasma parameters and forward scatter- 
ing of fast-electron bremsstrahlung. This section introduces 
these two models. The first is a fast-electron bremsstrahlung 
radiation model based on the Geant4 code, where the distribu- 
tion of X/7-ray sources refers to the shape of HL-3 plasma, 
and the X/y-ray detectors refer to the 2D HXR-imaging sys- 
tem. Other components, such as the stainless-steel vacuum 
chamber and tungsten copper shields, are also considered in 
this model. The second is the fast-electron forward-scattering 
model, in which the HXR emission strongly peaks in the for- 
ward direction, and the higher the electron energy, the nar- 
rower the emission cones. 

The simulation flowchart is shown in Fig. 3. We used 
an HL-3 tokamak as the particle-source configuration and 
an HXR-imaging system as the energy-deposition configura- 
tion. The particle-source distribution refers to the experimen- 
tal data from HL-3 discharges, and the particle-motion direc- 
tions are constructed using a fast-electron forward-scattering 
model. A sufficiently large number of photons were inputted 
into Geant4 and the code was run. 
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Fig. 3. (Color online) Simulation flow of X-ray photons energy de- 
position on the detector. 


Fig. 4. (Color online) Half-cutaway view of HL-3 device for Geant4 
simulation, where the green lines are X-rays, the brown stripes are 
toroidal and poloidal field coils, and the blue shell is the boundary 
of HL-3 plasmas. 


A. Geant4 Simulation Model 


We developed a bremsstrahlung simulation model based on 
the Geant4 code [42, 43] to simulate energy deposition on 
the detector. The geometry is defined as a considerably sim- 
pler structure than that of a realistic HL-3 tokamak, as shown 
in Fig. 4, which includes only a copper magnetic-field coil 
system and stainless vacuum vessel. An LYSO crystal and a 
single-energy X /y-ray are used in our simulations. The X /+7- 
ray source profile is assumed to have a Gaussian distribution 
in both the radial and vertical directions and has the same di- 
mensions as HL-3 plasma with R = 1.78 m, a = 0.65 m, 
and « = 1.3. The directions of all the electron motions are 
parallel to the magnetic axis. Geant4 tracks all X/y-ray pho- 
tons until they are absorbed in the materials and records the 
deposit energy and location when the X/y-ray photon enters 


the LYSO detector. 


B. Fast-Electron Bremsstrahlung Model 


The simulation is based on the forward scattering of X/7- 
ray bremsstrahlung for fast electrons [44], where all electrons 
are assumed to move from left to right, and each electron 
emits X/-ray photons with energies smaller than their ki- 
netic energies. The calculated results are shown in Fig. 5, 
where the X/y-ray emission cones are plotted for electrons 
with energies of 20, 50, 100, 200, and 500 keV. As shown 
in this figure, the higher the electron energy, the narrower 
the bremsstrahlung cones. Therefore, in a tokamak device, 
if the energetic electrons exhibit speed primarily along the 
magnetic-field lines, the X /7-ray emission is strongly peaked 
in the direction of the toroidal field and opposite to the direc- 
tion of the plasma current. 
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Fig. 5. (Color online) Bremsstrahlung-emission cones for mildly 
relativistic electrons. 


The forward-scattering characteristics are shown more in- 
tuitively in Fig. 6. Electrons (yellow dots) moving along the 
magnetic-field line (red dashed line) emit photons in all di- 
rections (green arrows), whereas the photon distribution is the 
densest (blue dots) in front of the motion. Fig. 7 can be used 
to describe the entire process: when the electron (red dot) 
moves to point A along orbit AA’, it has the movement di- 
rection of AB, the bremsstrahlung intensity is the strongest 


along AB, and the distribution characteristics are F (u, v, w) 
in all directions. 

The first step is to construct a particle-distribution func- 
tion. By adding three small variables 6,1,652,andé653 to 
each coordinate of A(Xo, Yo, Zo), we obtain a new point 
As(Xo +651, Yo +852, Zo +053). If 651, 652, andd,3 take dif- 
ferent values and form the set of {{051}, {6s2}, {6s3}}, and 
they are added to Xo, Yo, Zo respectively, we obtain a point 
set {{Xo+s1}, {Yotds2}, {Zo+653}}, which is centered on 
(Xo, Yo, Zo), and has a length of maz(és1) —min(d¢1) in the 
X direction, width of max(6,2)—min(0d.2) in the Y direction, 
and height of maz(6,3) —min(6.3) in the Z direction. There- 
fore, the main task in building a particle-distribution model is 
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Fig. 6. (Color online) Top view of the fast-electron forward- 


scattering diagram, where the red shaded area is the HL-3 plasma, 
the red dashed line is electron motion trajectory, the green arrows are 
the bremsstrahlung direction, and the blue dots are bremsstrahlung 
photons. 
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Fig. 7. (Color online) Fast-electron forward-scattering model, where 
point O is the center of HL-3, R is the minor radius of HL-3, point 
A is the location of electron, point B is in the movement direction 
of the electron, F (u, v, w) is the radiation direction function. 


to determine suitable values for 6.1, s2, andd.3. 

Referring to the HL-3 experiments with RF waves, the 
electron temperatures and densities approximately obeyed 
normal distributions in the radial and polar directions (left 
panel in Fig. 8). Therefore, only a set of random numbers 
651,052, andôsz that obey a normal distribution need to be 
constructed to obtain a distribution model similar to the elec- 
tron temperature or density profiles (right panel in Fig. 8). 

Notably, the direction of electron motion is parallel to 
the magnetic axis, that is, perpendicular to the cross section 
where the electrons are located. The electron-motion orbit is 
on an arc with radius R, such that one set of random numbers 
with a normal distribution can generate a particle distribution 
in the X and Y directions. Moreover, because the random 
distribution of particles in the Z direction is independent of 
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Fig. 8. (Color online) Left, the temperature data from HL-3 dis- 
charges (the blue dots and red stars are measured data, and blue 
curve and red curve are fitted profiles); right, the normal distribu- 
tion of incident particles on the HL-3 plasma cross section in Geant4 
simulations. 
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Fig. 9. (Color online)Schematic of fast-electron bremsstrahlung an- 
gle. 


the X and Y directions, the distribution of particles in the Z 
direction can be defined by another set of random numbers 
with a normal distribution. Additionally, the particle moves 
along orbital AA’, and the probability of its appearance at 
any point in the orbit is random, which can be achieved by a 
random angle 0. Let 


0 = irl ay (3) 


where £4 is a random number and 0 < €4 < 1. 

Two random numbers €,; and €,2 can be used to obtain 
a pair of random variables that follow the standard normal 
distributions N, and Ms, using the following formula: 


Ms = (21n€51) ` cos(27&.2) 
; (4) 
Ns = y (2ln€s1) « sin(27€ 52) 
where 0 < €51,€52 < 1. 
According to previous descriptions, let 
dst = 552 = Ms 
, (5) 
3 = Ns 


such that the point set {{Xo+6s1}, {Yo +852}, {Zo+6s3}} 
follows the standard normal distribution, and can be written 
as follows: 


Xs = 051-551 + Xo 
(051 -s1 + R) - cos(0) 
(—2ln€.1) - cos(27&,2) + R) - cos(@) 


= (O51 z 


Ys = 052 - 052 + Yo 
= (Ts2 - O52 + R) : sin(@) 
(—2ln€.1) - cos(27€s2) + R) - sin(@) 


, (6) 


= (752 ' 


Zs = 053° K: 653 + Zo 


= 053° K+ V/(—2lnE51) « sin(2rés2) + Zo 


where 051, 052, and gs3 are the plasma shape parameters 
and « is the elongation of the HL-3 tokamak. The particle- 
distribution function can be expressed as follows: 


G = G(Xs, Ys, Zs). (7) 


The next step is to construct an emission vector for each 
particle. When the fast electron moves along the orbit (red 


line) to A, its direction of motion is AB because the center of 
the electron motion is perpendicular to the Z axis; thus, A 
1 OŻ 
the direction of AB, and the probability of other directions 
follows the distribution function F(u, v, w). 

Assuming IJAB] = 1, the coordinate value of B(X,Y, Z) 
can be written as follows: 


. The radiation of fast electron A is strongest along 


X=X0o-Ý/R 
Y = Yy + Xo/R. (8) 


Z=% 


To simplify our simulation model, we add an auxiliary line 
CD to describe the process, where C.D passes through point 
B and is parallel to OA, as shown in Fig. 7. Similar to the pre- 
vious description of the particle-distribution model, if we add 
a set of small variables {{6a1}, {da2}, and{da3}} that follow 
the normal distribution to each direction of point B, and then 
project all the points onto the line of CD, we obtain a point 
set {{X + da}, {Y + 6a2},{Z + 6a3}}, which follows the 
normal distribution along the line of CD and is located in 
the region of A, A», as shown in Fig. 9, where the maximum 
emission angle is P. The relationship between the maximum 
emission angle and length of A, Ag is 


—> 
AA 
tan() = = (9) 


The emission angle is obtained by varying the lengths of 
Aj, Ao. 


Point A is connected to each point on the set {{X + 
dai}, {Y + ba2},{Z + 6a3}}, and these lines are considered 
as the propagation direction of the X/ y-ray photons. Thus, 
we obtain a forward-scattering model for each electron. 

Two random numbers €4; and €49 are used to obtain a pair 
of random variables M4 and Ng, which follow a standard nor- 
mal distribution: 


Ma = y (2ln€a1) - cos(27€a2) 
, (10) 
Na = y (2ln€a1) - sin(27€a2) 
and then let 
ba1 = a2 = Ma 
(11) 
6a3 = Na 


Point set {{X + dai}, {Y + da2},{Z + 6a3}} follows a 
standard normal distribution and can be written as follows: 


Xp=X+ (oa - dai) ` Yo/R 
= X + V/(—2ln€a1) - cos(2m€a2) - Yo/R 


Yp=Y+ (oaz - Oa2) -Xo/R 


. a2 
= Y + V/(—2ln€a1) - cos(2ma2) - Xo/R i 


Zp = Z + 043° ba3 
= Z + /(—2ln€a1) - sin(27€a2) 


where 0 < €a1,€a2 < 1 and oai, Oa2, Cag are the emission- 
angle parameters. The emission vectors of each particle are 


u = Xp — Xo 
= y (—2ln€a1) - cos(27Ea2) - Yo/R 


v= Yp — Yo 


; (13) 
= y (—2lnéa41) - cos(2mEa2) + Xo/R 


W = Zp ia Zo 
= y (—2lnéa1) - sin(27€a2) 
The emission-angle distribution function can be expressed 
as 
F = F (u,v, w). (14) 


Because the cross section of the X-ray emission cone is 
circular, that is, symmetrical in the X, Y, and Z directions, we 


obtain oq = Ca2 = Caz. If the emission angle ® is known, 


we can derive A; A> from Eq. 9 and the values of c41, Caz, 
and oq3. 

In summary, the key to the construction of the fast-electron 
bremsstrahlung model is to determine the particle-distribution 
and radiation-direction models. The particle-distribution 
model G(Xs, Ys, Zs) generates the location of each photon, 
and the radiation-direction model F (u, v, w) indicates the di- 
rection in which the photons move. 


IV. SIMULATION RESULTS 


Geant4 traces each X /y-ray photon until it is deposited on 
the target detector or lost from space, and the code records all 
particle-deposition locations and deposition energies on the 
detector. The output data format is listed in Tab. 1. 


Fig. 10. (Color online) Simple schematic of the detector, where the 
blue squares are detector units, the X-, Y-, and Z-axes represent the 
thickness, width, and height of the detector, respectively, and dy and 
dz are the distance of every two adjacent detector units in the Y and 
Z directions. 


TABLE 1. Data structure of Geant4 simulation outputs. 


X Y Z Deposition Energy 
zı yı zı Eı 
T2 y2 z2 E2 


X,Y, andZ are the thickness, width, and height of the de- 
tector, respectively. (As shown in Fig. 10, each blue square 
is a detector pixel, and the distance between each pair of ad- 


jacent detectors are d, and d, in the Y and Z directions, re- 
spectively. 

Because the size of the deposition position of the individual 
ray photons is almost negligible, the resultant radiation image 
consists of independent discontinuous pulses. To reconstruct 
the radiation image, the detector must be separated into uni- 
form intervals and the deposited energy in each interval must 
be integrated. The main part of the algorithm is as follows: 


s(i, j) = sum(E( y > Ymin + (i — 1) x dy & 
Y < Ymin +1 X dy & 
z > Zmin + (j — 1) x dz & 
z < Zmin + j X dz)) 


(5) 


where s is the deposited energy in each unit area, ¿ and 
j represent the cyclic variables in the Y and Z directions, 
Ymin and Zmin represent the minimum values of Y and Z, 
respectively, and d, and d, are the step widths in the Y and 
Z directions, respectively. 


A. Cases with Different Particle Numbers 


We first consider the same size of the plasma cross section 
(plasma radius of 0.65 m) and gradually increase the number 
of incident photons. Let osı = 052 = 053 = 1/6. We can 
obtain the particle distribution in the X,Y, and Z directions 
using the particle-distribution model from Eq. 6 as shown in 
(a), (b), and (c) in Fig. 11. The emission angle ® = 150° and 


length of A; Aə can be calculated using Eq. 9; therefore, we 
can obtain o@1, 0d2, a3 and generate the direction distribu- 
tion along CD from Eq. 12, as shown in Fig. 11(d). Fig. 12 
provides a more intuitive distribution map of the particle dis- 
tribution in the HL-3 tokamak, where the color represents the 
normalized particle density, with the warmer colors indicating 
denser particles and cooler colors indicating fewer particles. 
The outer dark grids represent the plasma boundaries. 

The simulation results are shown in Fig. 13. In this fig- 
ure, the color bar represents the normalized energy-deposition 
value: the warmer colors represent more energy deposition, 
and the cooler colors represent less energy deposition. (a), 
(b), (c), (d), (e), and (f) represent the incident photon numbers 
of 2x 108, 4x 108, 1x 10°, 5x 109, 1x 10", and 2x 107°, re- 
spectively. When the number of incident photons is small, the 
energy deposition is relatively dispersed, and no obvious dis- 
tribution characteristics are observed. However, as the num- 
ber of incident photons increases, the energy-deposition dis- 
tribution gradually becomes concentrated in the middle area, 
similar to the particle distribution shown in Fig. 12. 


B. Cases with Different Plasma Sizes 


We consider different sizes of plasma cross sections, all 
with sufficient incident photons of 2 x 101° and fixed emis- 
sion angles of 150°. The particle-source size can be adjusted 
by changing the values of 051,052, and 0,3. Fig. 14 shows 
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Fig. 11. (Color online) Normalized particle distribution (a) in the 


X direction, (b) in the Y direction, (c) in the Z direction, and (d) 
particle distribution around point B at ® = 150°. 
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Fig. 12. (Color online) Normalized particle distribution on plasma 
cross-section, where the warmer colors indicate denser particles and 
cooler colors indicate fewer particles; the outer dark grids are plasma 
boundaries. 


the cross section of the HL-3 tokamak and plasmas, in which 
(a), (b), (c), (d), (e), and (f) exhibit particle-source radii of 
0.1m, 0.2m, 0.3m, 0.4m, 0.5m, and 0.65m (the maxi- 
mum size of the HL-3 plasma), respectively, where the cor- 
responding plasma-size parameters (051 = Os2 = 053) are 
1/40, 1/20, 1/10, 1/5,1/2, and 1/1, respectively. The color 
in each figure represents the normalized particle distributions, 
where warmer colors indicate denser particles and cooler col- 
ors represent fewer particles. 


Fig. 15 shows the simulation results. Even though the 
energy-deposition areas vary with different plasma cross- 
section sizes, all are symmetrically distributed. This suggests 
that the energy-deposition characteristics on the detectors are 
independent of the plasma size. 


X (mm) 


Fig. 13. (Color online) Normalized energy deposition on the detec- 
tors, where the warmer colors indicate higher energy deposition and 
cooler colors indicate lower energy deposition. (a) to (f) are incident 
photon numbers of 2 x 10°, 4 x 108, 1 x 10°, 5 x 10°, 1 x 107°, 
and 2 x 10'°, respectively. 


C. Cases with Different Particle-Emission Angles 


In HL-3 tokamak discharges, the energy of the ray pho- 
tons generated by fast-electron bremsstrahlung is continu- 
ous. However, in our simulations, all initial photons are mo- 
noenergetic. According to the description in Section III, the 
cone angles of photon emissions generated by fast-electron 
bremsstrahlung with different energies are different; thus, we 
can use different cone angles to represent electrons with dif- 
ferent energies. 

The following formula assist in understanding this process: 


N 
E=k v’ (16) 
where F is the total deposited energy of the incident photons 
on the detector, k is the energy of each photon, JN is the total 
number of photons, and V is the space volume of the photon 
emission. 

Because k is a constant value in our simulation, if the to- 
tal incident photon number NV is the same for each calcula- 
tion, then the total energy E can be adjusted by changing the 
photon-emission volume V. Eq. 16 shows that the smaller 
the emission volume V, the larger the value of N/V, and the 
higher the energy deposition, and vice versa. The radiation 
model in Fig. 5 shows that the higher the electron energy, 
the narrower the bremsstrahlung cones, and as a result, the 
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Fig. 14. (Color online) Normalized particle distribution on plasma 
cross section, where warmer colors indicate denser particles and 
cooler colors indicate fewer particles. (a) to (f) exhibit plasma cross 
sections with sizes of 0.1, 0.2, 0.3, 0.4, 0.5, and 0.65 m, respectively. 


smaller the emission volume. Hence, the electron energies 
are related to the radiation-cone angles; a larger cone angle 
represents a smaller electron energy, whereas a smaller cone 
angle represents a larger electron energy. 


In this subsection, we investigate the influence of the 
particle-radiation angle on energy deposition, all cases have 
sufficient incident photons of 2 x 101° and a fixed plasma ra- 
dius of 0.65m. Fig. 16 exhibits four different radiation angles. 
To display all the angles in the same diagram, we divide the 
cone radius of the larger angle by the different coefficients. 
The red cone angle represents a radiation angle of 60°, twice 
the green cone radius represents a radiation angle of 120°, 
three times the blue cone radius represents a radiation angle 
of 150°, and six times the black cone radius represents a ra- 
diation angle of 170°. 

Fig. 17 shows the simulation results under these radiation 
angles, (a), (b), (c), and (d) represent radiation angles of 
60°, 120°, 150°, and 170°, respectively. The color bar in 
the figure represents the normalized energy-deposition value, 
where warmer colors represent more energy deposition and 
cooler colors indicate less energy deposition. A white ellip- 
tical dashed line is used to label the shape of the energy de- 
position. The diagram shows that when the radiation angle 
is small, the shape of the energy deposition tends to be cir- 
cular, and when the radiation angle is large, the shape of the 
energy deposition becomes elliptical, and the aspect ratio of 
the ellipse increases with the angle. However, in all cases, the 
deposited energy is symmetrically distributed. 

Although single-photon energy is used in our calculation, 
changing the radiation angle can represent different energy 
values to some extent. From the calculation results, changing 
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Fig. 15. (Color online) Normalized energy deposition on the de- 
tectors, where warmer colors exhibit higher energy deposition and 
cooler colors indicate lower energy deposition. (a) to (f) show cases 
of different plasma cross sections from Fig. 13. 
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Fig. 16. (Color online) Fast-electron emission angles, where the red 
lines represent a maximum emission angle of 60°, the green lines 
represent a maximum emission angle of 120°, the blue lines repre- 
sent a maximum emission angle of 150°, and the dark lines represent 
a maximum emission angle of 170° 


the radiation angle of the photons does not alter the symmetry 
distribution of the energy deposition. We believe that similar 
results can be obtained if the energy of the particle and its 
associated radiation angle are used as variables. 
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Fig. 17. (Color online) Normalized energy deposition on the detec- 
tors, where the warmer color represent higher energy deposition and 
cooler colors indicate lower energy deposition. (a) to (d) show cases 
of different emission angles from Fig. 16. 


V. ANALYSIS MODEL OF ENERGY DEPOSITION ON 
THE DETECTOR 


A 2D analytical model is developed to further explain the 
energy deposition that appears on the detector. Two main rea- 
sons exist for the symmetrical distribution of the deposited 
energy on the detector: the first is the forward-scattering ef- 
fect of bremsstrahlung, and the second is the pinhole imaging 
effect. Owing to the symmetry in the Z direction, we con- 
sider only the middle plane of the HL-3 plasma. In Fig. 18, 
the red region is the midplane plasma, the black arcs are the 
plasma boundaries, and the green rectangle on the left is the 
location of the HXR-imaging system. 


Fig. 18. (Color online) Middle plane of HL-3 plasma (red shadow), 
where the red dashed line is the magnetic axis and the left green box 
is the HXR-imaging system. 


Energy deposition on the detector from a 2D plasma plane 
can be calculated by separating the plasma surface into mul- 
tiple lines. For convenience of calculation, we list only a few 
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Fig. 19. (Color online) Normalized energy-deposition probability 
on the detector. Left, particles from HOPO, H1P1, H2P2, H3P3, and 
H4P4; right, average energy-deposition probability from these five 
lines. 


lines, as shown in Fig. 18 with HOPO, H1P1, H2P2, H3P3, 
and H4P4; all the lines pass through the center point, and the 
angle between each pair of adjacent lines is 30°. The number 
of particles between H* and P> in each line follows a Gaus- 
sian distribution (consistent with the radial distribution of the 
electron temperature or density). 

As introduced in Sec. II B, the probability of the photons 
emitted by an electron bremsstrahlung follows a Gaussian 
distribution centered on its direction of motion. The collima- 
tion pinhole in an HXR-imaging system is small and regarded 
as a point. Based on the principle of pinhole imaging, only 
photons in the direction of the connecting line between the 
positions of the electron and pinhole can reach the detector. 
If we record the deposited energy and location of each pho- 
ton and add them, we obtain the energy distribution on the 
detector. Taking the i-th electron as an example, the position 
of the energy deposition can be calculated using a triangular 
relationship, and the energy value of the deposition can be 
calculated using the following equation: 


E; = Ep X Pai X Pei, (17) 


where E; is the normalized energy, Æp is the energy of a 
single photon, pq; is the normalized probability of electrons 
appearing at that position, and pe; is the normalized probabil- 
ity of a photon passing through the pinhole. 

Because the particle source is linear, the energy deposited 
on the detector is also linear. The calculated results are shown 
in Fig. 19. The deposited energy of photons emitted by 
electrons on the HOPO line on the detector is symmetrically 
distributed about the middle line (the black curve in the left 
figure), whereas the deposited energy from other lines devi- 
ates from the middle line (the green, blue, red, and purple 
lines in the left figure). However, as shown in the diagram, 
the deposited energy from H1P1 and H3P3 are symmetri- 
cally distributed about the middle line, which is the same 
energy-distribution characteristic from lines H2P2 and H4P4. 
Adding the calculation results from all lines shows that the 
total deposited energy on the detector is still symmetrically 
distributed (the right figure). 

Although our analytical model is 2D, we believe that if we 
consider the entire plasma range, a symmetric distribution of 
energy deposition still exists. 


VI. SUMMARY AND DISCUSSION 


A fast-electron bremsstrahlung imaging model is devel- 
oped using the Monte Carlo simulation software Geant4, 
which is based on the plasma geometry and HXR-imaging 
system at the HL-3 tokamak. The distribution of fast elec- 
trons refers to the temperature and density during the RF wave 
cutrent-drive experiments, and the forward scattering of fast- 
electron bremsstrahlung radiation is considered. Preliminary 
simulation results show that the energy deposition on the de- 
tector is symmetrically distributed as long as the number of 
particles is sufficiently large, even if the toroidal effect causes 
an asymmetrical plasma distribution. The deposited energy 
on the detector appears to have a distribution similar to that 
of the electron temperatures or densities of the plasma cross 
section during the RF-wave experiments. Multiple simula- 
tion results confirmed that the symmetric distribution char- 
acteristics remained when the particle-distribution model or 
fast-electron bremsstrahlung angle was changed. 

Two-dimensional HXR-imaging systems are particularly 
important for the study of fast-electron physics; however, the 
spatial asymmetry of plasma poses some difficulties in an- 
alyzing experimental data. The simulation results compen- 
sated for this deficiency by constructing a relationship be- 
tween the energy deposition on the detector and spatial dis- 
tribution of the plasma parameters, indicating that the diag- 
nostic data from the 2D imaging system can be directly used 
to analyze the 2D spatiotemporal evolution properties of fast 
electrons. Moreover, if the energy deposition on the detector 
is known, we can deduce the spatial-distribution character- 
istics of fast electrons or the fast-electron forward-scattering 
properties. 

For future fusion-reactor-scale devices, the higher the 
plasma parameters, the greater is the role of fast electrons 
in maintaining the steady-state operation of the system. Ac- 
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